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Abstract. We determine projected rotation velocities v sin i in DAZ white dwarfs, for the first time using the 
rotational broadening of the Call K line. The results confirm previous findings that white dwarfs are very slow 
rotators, and set even more stringent upper limits of typically less than 10 km/s. The few exceptions include 3 
stars known or suspected to be variable ZZ Ceti stars, where the line broadening is very likely not due to rotation. 
The results demonstrate that the angular momentum of the core cannot be preserved completely between main 
sequence and final stage. 
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1. Introduction 

The last years have seen a new generation of stellar evolu- 
tion models emerging, which try to incorporate the influ- 
enc e of rotation on internal structure and evolution (see 
e.g. iMev net fe MaederLll997HLanger et ailll998l:lLangeil 
ll998HHeger et al.ll2000l and later papers of these groups) . 
These studies have been motivated by the realization 
that stellar surface abundances in many areas of the 
Hcrtzsprung-Russell diagram require mixing processes in 
addition to convection. The effects o f rotation in mas- 
sive s tars are described in detail in iMaeder fe Mevnetl 
(2000j), which gives also many relevant references. In the 
calculation of rotating stellar models inevitably the ques- 
tion of angular momentum transport arises, which can- 
not yet be described self-consistently from first principles. 
Ordinary molecular viscosity is much too small, and other 
mechanisms have to be in voked, e.g. turbu l ence c reated 
through sh e ar in stability ilSpiegel fc ZahnL Il970t IZahnL 



I1974L 1197.1 |l992t) generation of internal gravity w aves 
l|Talon et all 12002 : iTalon fc Charbormel I2003L 1200^ . or 



magnetic torques l)Spruitl ll998h 

In this situation the observation of white dwarf ro- 
tation can provide some constraints since it tells us the 
remaining angular momentum of the typical 0.6 M inte- 
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rior core of the progenitor at the end of stellar evolution. If 
the total angular momentum of the core is preserved one 
would expect equatorial velocities in the range of 15 km/s 
for solar type slow-rotating progenitors to 2000 km/s for 
a fast-rotating AO star, corresponding to rotation periods 
from 1 h down to a few minutes (a number useful to re- 
member is that for a 0.6 M Q white dwarf an equatorial 
velocity of 15 km/s corresponds to a 1 h rotation period). 

The only detailed calculation with a sophisticated 
treatment of angular momentum transfer, which makes a 
specific prediction f or the fi nal whi te dw arf rotational ve- 
locity is the work of lLanger et~aT] l)l999|) . For the 0.8 M Q 
white dwarf resulting from the evolution of a rotating 
3 M Q model they expect an equatorial rotation of 28 km/s. 

On the other hand - assuming the core to remain 
strongly coupled to the outer expanding envelope in the gi- 
ant stage - the rotation of the final white dwarf should be 
extremely slow , with periods in the range 30 to 1000 years 
l|Spruid . Il998|) . Rotation periods of hours to days have in 
a few cases been determined from the rota tional splittin g 
of pulsation frequencies; a list is given in ISpruitl 1^)98), 
to which could be a dded a p e riod o f 13 hrs for L19-2 and 
58 hrs for GD165 l)BradlevL 12001)) . Periods around one 
day have also been measured from variable pola r izatio n 
ilSchmidt fc NorsworthvL Il99ll ISchmidt fc Smith! . Il995lh 
but for some of these magnetic white dwarfs no variation 
has been observed, implying rotation periods of at least a 
century. 
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Fig. 1. Call K line profile in WD 1826-045. Dotted: ob- 
served profile; continuous: non-rotating model after con- 
volution with instrumental profile (falls on top of the ob- 
servation in the line center!); dashed and dot dashed with 
decreasing line depths: rotating models with v sin i = 5, 
10, 20 km/s 



The great majority of rotation data for white 
dwarfs have been determined spectroscopically, exclu- 
sively using the sh arp NLTE core of th e Ha li ne in 
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The sample now contains about 50 stars with determined 
vs'mi. A few of these objects show a line broadening sug- 
gesting velocities of a few dozen km/s, but in many cases 
the unsatisfactory line fits indicate that the reason is not 
rotation but magnet ic splitting or unknow n, as in the case 
of the ZZ Ceti stars IjKoester et allll998j) . In the majority 
of objects only upper limits can be determined which typ- 
ically range from 10 to 30 k m/s and are thus ma rginally 
compatible with the result of lLanger et alJ 1)19991) . While 
this does not rule out t hat mos t whit e dwarfs could be 
extremely slow rotators, pprui 3 <|1 998ft . taking all avail- 
able empirical evidence and following his idea of strong 
core-envelope coupling felt a need for mechanisms, which 
would accelerate the white dwarfs again to the observed 
values after the efficient loss of angular momentum in the 
giant phase. 

There is thus clearly a strong motivation to deter- 
mine more stringent upper limits using the spectroscopic 
method. This has become possible after the detection of 
a large number of DAZ objects found in th e stud ies of 
IZuckerman fc Reidl l(l998h : IZuckerman et alJ l|2003h and 
as a byproduct of t he search for Supernova la proge nitors 
l|Napiwotzki et all I2003L l200lHKoester et alll2005ft . The 
Call K line in these hydrogen-rich white dwarfs is deep 
and quite narrow and should be useful for measuring much 
lower rotation velocities than is possible with Ha. In this 
paper we report the results of such a study. 



velocities of white dwarfs 

Table 1. Atmospheric parameters T c g in K, surface grav- 
ity log g, and Ca abundance for the sample. [Ca/H] is 
the logarithm of the Ca/H abundance ratio by numbers. 
C olumn Ref gives the sou rce o f the parameters wi th Z 
= IZuckerman et all fcOO.Sh . R = lRollenhag"en~l l(2004l) and 
IKoester et all l|2005|) . Columns K and S give the number 
of available spectra from the Keck or the SPY samples. 
Please note that the Ca abundance for HS 0047+1903 in 
\Koester et all \200A) was based on an erroneous equiva- 
lent width. The correct abundance is given in Tabled 
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2. Observations 

The sample consisted of spectra for 38 single cool DA 
white dwarfs coming from two different sources. The first 
is the SPY survey (ESO Supernova la Progenitor Survey) 
for double-degenerates. In the search for variable radial 
velocities a large number of high-resolution spectra are 
obtained with the UVES spectrograph at the ESO VLT 
Paranal Observatory, which ar e usef ul also for man y other 
projects. iRollenhagen I l)2004) and IKoester et al.1 l(2005h 
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used this sample to identify hydrogen-rich white dwarfs 
with Call lines and those objects form the first part of 
our present sample (excluding those where the lines are 
interstellar). Details of the obs ervations and the reduc- 
tion procedures ar e descr ibed in iNapiwotzki et all l)2003(l 
and iKoester et all l|200lj) : we therefore repeat here only 
the facts relevant for this project. The observations used 
a wide slit of 2'/l, because the SPY was conducted as a 
filler project for bad weather conditions. This is important 
since the nominal resolution with this slit is 18500 corre- 
sponding to 0.21 A at the Call K line, however, the actual 
resolution can be higher if the seeing disk is smaller than 
the slit. The actual seeing during the SPY observations 
typically varied between l'.'O to l'/5. 

The second source of sample objects are Keck obser- 
vations of DAs in a search f or DAZs with meta l s. The 
obs ervations are described in IZuckerman fc Reidl Jl998) 
and lZuckerman et al.l l|2003j) . Here the nominal resolution 
with a slit width of 1'.'15 is 34000, corresponding to 0.12 A 
at Call K, but also in this case the actual resolution might 
be slightly higher at times of excellent seeing. 

For most objects more than one spectrum is available, 
for some even spectra from both surveys. Atmospheric pa- 
rameter s T P ff , log q, and Ca a bun dances for the sample are 
given in lKoester et all l)2005l) and lZuckerman et all ((2003) 
(which also show some representative spectra) and are col- 
lected in Tabic n For our analysis the spectra were shifted 
slightly to bring the K line in agreement with the labora- 
tory wavelength and then normalized to a continuum of 
1.0 with a straight line fitted to the continuum approxi- 
mately 5 A on both sides of the line center. 



3. Determination of rotation velocities 

were determined 



theoretical 



Rotation velocities were determined using 
model atmospheres with spectral lines broadened by a 
rotational broadening function (see e.g. lUnsoldl ^)68). 
Our model grid is calculated assuming LTE (local ther- 
modynamic equil ibrium) with methods and input physic s 
des cribed inlFinley et alJ l)l997l) . lKoester fc Wolfj l|2000|) . 
and lHomeier et alJ l|l998f) . The Call lines have the advan- 
tage compared to the Balmer lines - the only lines used in 
the past for spectroscopic determinations in white dwarfs 
- that the central Doppler core from thermal broadening 
is much narrower because of the higher atomic weight. 

On the other hand, the Balmer lines originate from 
the dominant element in DAs and therefore only two pa- 
rameters, T e g and log g, are needed to describe the fitting 
model. In the case of the K line an addit ional parameter is 
the C a ab undance. The results given in IZuckerman et al.1 
l)2003|) and lKoester et alJ l|2005|) have been determined us- 
ing interpolations in pre-calculated tables of equivalent 
widths. In the calculations of the theoretical equivalent 
widths a fixed interval was assumed, and no attempt was 
made to match this exactly in the observed spectra; more- 
over, the spacing of the grid of abundances was rather 
coarse. It is therefore not surprising that calculated spec- 
tra for the "exact" atmospheric parameters using these Ca 
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Fig. 2. Typical examples for the determination of rota- 
tional velocities for the Call K line, ranging from some 
of the best spectra to marginal. Dotted: observed profile; 
dashed: non-rotating model after convolution with instru- 
mental profile; continuous: best fit with rotation. In the 
two top figures the best fit is zero rotation and the models 
fall on top of each other. 

abundances did not always match the spectra perfectly. 
We have then changed the Ca abundance (typically less 
than a factor of two) until the EW of the observation 
agreed exactly with that of the theoretical model. This 
is obviously a prerequisite for the determination of rota- 
tional broadening of the profile. 

The standard rotational broadening profile uses two 
free parameters, the limb darkening coefficient (3 and the 
projected equatorial rotation velocity v sin i. Here i is the 
inclination of the rotation axis against the line-of-sight, 
ranging from to 90 degrees (note that only this combi- 
nation can be determined with the spectroscopic method) . 
Our theoretical models show that f3 varies from w 0.2 at 
the center of the K line to « 1.2 in the far wing. Since the 
parameter has almost no influence on the result we have 
continued to use (3 = 0.15 which is used by all authors for 
Ha. The best fit value for v sin i and the errors were then 
determ ined with a y 2 fittin g routine very similar to that 
used in lKoester et alJ l)l998|) . Figure ^ shows the observed 
line in WD 1826-045 as a typical example together with 
theoretical models with different rotational velocities ap- 
plied. The figure demonstrates that rotational velocities 
of a few km/s can easily be measured from the Ca line. 

3.1. A remark on the standard spectroscopic method 

Assuming that the intensity of the radiation field at the 
surface of a star depends only on wavelength A and the 
angle 9 against the normal, the observational quantity for 
a non-rotating star is the disk-averaged intensity 



J.I 

7T 



A) dxdy 



where [i = cos 8, and x, y are Cartesian coordinates on the 
visible surface of the star in units of the stellar radius R. 
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In the case of a rotating star we have in addition a 
Doppler shift, which depends on the position on the disk 

+l +Vi-x' 2 

J=i J J A\(x,y))dxdy 

If the y-axis is chosen parallel to the rotation axis, AA 
will depend on x only. This integral can easily be evalu- 
ated numerically if theoretical models are available for the 
intensity as a function of /i. It is common practice, how- 
ever, to proceed with an analytic evaluation by assuming 
a limb-darkening law 

7( M ) = /c((l-e)+eM) 

with central intensity I c . This law can also be written 
relative to the intensity at the limb J; 

!(/*) = 1,(1 + /?/*) 

with the limb-darkening coefficient (3 = e/(l — e) used 
above. Introducing this into the integration over the disk, 
with a function of position x, y, the integration can be 
transformed into a convolution of the disk-averaged inten- 
sity for a non-rotating star 

+i 

I rot = J G(x)I(\ - xAXi) dx 

-i 

with the maximum shift at the limb AA; , and 
n( . _ 4(1 - e) + eVl - x 2 + 7re(l - x 2 ) 

G[x) Mi - 6/3) 

fseee.g. lljnloldLll968h . This is a large simplification, since 
only the theoretical disk-averaged intensity (also called 
the astrophysical flux) needs to be calculated, which is al- 
ways needed for the analysis of spectra, rotating or not. 
However, the derivation makes the assumption that the 
limb darkening law is independent of wavelength, in par- 
ticular that the coefficients e or (3 are constant over the line 
profile. This is clearly not the case for the Call K line as 
mentioned above and it is also not true for Ha. We have 
therefore studied the effect of this assumption by com- 
paring the exact direct integration with angle-dependent 
intensities with the simplified method using the convolu- 
tion technique. While this has been studied before in the 
literature for other types of stars it has to our knowledge 
never been tested for white dwarf spectr a. 

The result of this test is llPruinl 120051 that small differ- 
ences between the two methods exist in the line cores, but 
that the error made when using the simple convolution 
method is much smaller than 1 km/s, and thus negligible 
in our study. 

3.2. Observed and theoretical line profiles 

When comparing the observed and theoretical profiles af- 
ter taking into account instrumental broadening with the 



nominal resolution in the latter, we obtained nearly per- 
fect fits for about half of the objects without rotational 
broadening, implying zero or at least extremely low rota- 
tion. In the other half of objects the observed line profile 
had a deeper core than our "unrotated" models. Since ro- 
tational broadening can only decrease the central depths, 
we could thus not obtain a fit. 

What could be the reason for this discrepancy? The 
line profiles are of course influenced by the reduction pro- 
cess, e.g. the subtraction of sky background, to mention 
just one step which could affect especially the line cores. 
Indeed, when we have several spectra of the same object, 
there are sometimes differences in equivalent width and 
also line depth. Another possibility is on the modeling 
side: the largest discrepancies (in seven objects) occur for 
very deep lines, where the core is saturated. In these cases 
the line core originates in the outermost layers of the at- 
mosphere and it is quite possible that NLTE effects (not 
included in our models) will make the model lines deeper, 
as is well known for Ha. 

We believe, however, that in most objects the reason 
is that the actual resolution of the observation was bet- 
ter than nominal. This is supported by the finding that in 
several spectra of the same stars the equivalent width is 
constant but the line profile changes, with the line being 
deeper and narrower in some of the spectra. This strongly 
indicates a changing actual resolution for periods with bet- 
ter seeing. 

Unfortunately we do not know the real resolution for 
the observation. In order to be able to nevertheless extract 
some information on the rotation velocities we have there- 
fore increased the resolution (narrowed the instrumental 
profile) until our models achieved a good fit without ro- 
tation. We have then used our \ 2 routine to determine 
v sin i. Obviously the result must be compatible with zero, 
but we nevertheless obtain a reasonable upper limit from 
the errors provided by the fit. The results obtained for 
31 objects for which this procedure led to a good fit with 
the nominal or within a plausible range of resolutions (see 
below) are shown in Table El A number in the resolution 
column indicates the higher than nominal resolution used 
in the fit. 

There is of course some arbitrariness involved with 
this method, since the resolution could actually be even 
higher, and could also have been higher in the cases where 
our zero-rotation models had lines deep enough. In order 
to analyze the uncertainty we have assumed a maximum 
plausible resolution of 39500 and 65500, equivalent to see- 
ing of l'.'O and 0'.'6, respectively, for SPY and Keck spectra. 
The column v sin i(max) gives the results of the fitting 
procedure, which should be considered as very firm upper 
limits. 

Fig El presents four typical examples from spectra with 
high to rather low signal-to- noise ratios. 
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3.3. Error estimates 

We estimated the errors caused by the uncertainties of 
T e ff, log g, [Ca/H], (3, and the resolution. The summed er- 
rors from the first four quantities were of the same order of 
magnitude as the formal la-error of the x 2_ fitting given in 
Table |21 The greatest error was caused by the uncertainty 
of the resolution and can be estimated by comparing the 
values in the column v sin i and (v sin i) max . 

3.4. Objects with very deep line cores 

Even at the assumed maximum increase of the nominal 
resolution there remain seven objects which cannot be fit- 
ted by unrotated profiles. All of these objects have strong 
Ca lines, which are saturated in the core. In these cases 
we have used a Lorentz profile with zero central intensity 
and a width calculated to reproduce the observed EW. 
These profiles were then used to fit the observed profiles, 
resulting in firm upper limits to the rotational velocities. 
These results are given in Table The spectra and fits 
are shown in Fig. |3J where the dotted lines are the obser- 
vations, the dashed lines Lorentz profiles with zero central 
intensity as an extreme simulation of the line profile, and 
the continuous lines are the "rotated" profiles. 

4. Discussion 

21 objects in Table El and all 7 in Table are com- 
patible with zero rotation velocities (within ler errors). 
Even allowing for a higher than nominal resolution, within 
the plausible range, gives upper limits typically smaller 
than 10 km/s. This is a clear improvement compared 
to the spectroscopic determination from the core of Ha, 

10 - 30 km/s 
nfortunately we 
are aware of onl y two Ha determinati ons for objects of our 
sample, both in lKoester et~aTl l|l998h . For WD 2115—560 
the upper limit from Ha is 35 km/s, while we obtain 
11 km/s from the Ca line. The other common object is 
the ZZ Ceti star WD 2326+049 (G29-38). From Ha a ve- 
locity of 45 ±5 km/s is obtained, but the authors note that 
the fit is not good and the broadening very likely not due 
to rotation. From our Ca determination we find a range 
of 14-28 km/s from different spectra. 

G29-38 has been very extensively studied as a 
pulsator, both w it h ph otometry and spectroscopy, 
van Kerkwiik et"ai] {2000) detected periodic Doppler 
shifts in the line profiles with amplitudes up to 5 km/s, 
which they tentatively associ ated with the hori z ontal mo- 
tions of g-mode oscillations. iThompson et al. | il2003h re- 
cently found Doppler shifts in G29-38 corresponding to 
radial velocity changes up to 16.5 km/s. By combining a 
large number of spectra with close to zero Doppler shift 
and comparing with an average spectrum over several pe- 
riods they found no significant difference in the line pro- 
files. They concluded that pulsation cannot be the rea- 
son for the flat-bottom line profiles. This conclusion may 
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be somewhat premature, however, since the study only 
demonstrates that even the profile in one particular pul- 
sation phase does not fit a rotationally broadened model 
spectrum. It is plausible that even with an average shift of 
zero velocity in that phase different parts of the stellar sur- 
face may show positive as well as negative velocity along 
the line-of-sight, leading to zero shift but a broadening of 
the profile. This can only be decided with numerical sim- 
ulations of the pulsating surface including all shifts and 
realistic center-to-limb variations. In any case the differ- 
ence between the Ha and Ca results and the significant 
variation of the Ca "rotation velocities" confirm convinc- 
ingly that rotation is not the origin of the broadening. 

WD 2326+049 is one of the 10 objects for which we 
determine a non-zero rotation velocity for the best fit- 
ting model. Two others are WD 1150-153 and WD 1204- 
136, which according to their atmospheric parameters fall 
within or near the ZZ Ceti instability strip. WD 1204- 
136 was found to be non- variable (upper limit 0.3%) by 
iBergeron et~al1 l)2004h . while WD 1150-153 to our knowl- 
edge has not been observed for variability. The lat- 
ter is on our candidate list for upcoming time-resolved 
photometry, but in any case it seems very likely that 
the same mechanism - as yet unidentified - is respon- 
sible for the line s h apes as in the known ZZ Ceti 
stars i Koester et all Il998j) . Four other white dwarfs 
with non-zero rotation (WD 0543+579, WD 1407+425, 
WD 1344+106, WD 1821-131) are cool and the param- 
eters are based on photometry only. In the first two cases 
the parallax is not known and thus the s urface gravity can- 
not b e determined and was fixed at 8.0 (|Zuckerman et ail 
l2003h . However, even in these cases the rotation velocities 
are definitely very small, < 18 km/s at most. 

The results of this paper confirm the findings of all 
previous studies, that white dwarfs are slow rotators, and 
set even more stringent limits for those objects with Ca 
lines. Typical rotation velocities are below 10 km/s and 
may even be close to zero velocity for many objects. While 
this could be just marginally consistent with conservative 
evolution for solar-type slow-rotating progenitors, the re- 
sult would definitely exclude conservation of angular mo- 
mentum for more massive progenitors. How many of such 
white dwarfs could we expect in our sample? 

Assuming for simplicity a unique initial-final mass 
relation between r emnant and progenitor, we find from 
IWeidemannl Ij2000r ) that a 0.6 M Q white dwarf originates 
from a 2 M Q main sequence progenitor. 0.6 Mq (at T c ff = 
10000 K) correspond to \ogg = 8.0. Taking the numbers 
for the atmospheric parameters in Table^at face value, we 
have 27 objects with a logg determination (versus assum- 
ing 8.00), of which 11 have \ogg > 8.0 and should have 
originated from fast-rotating main sequence stars more 
massive than 2.0 Mq. While the exact numbers may be 
uncertain, there is no doubt that these results mean that 
angular momentum in the core cannot be pre served. These 
result s also clearly rule out the prediction bv lLaneer et alJ 
(1999J, which were based on stellar model calculations 
without magnetic fields. 
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While the original study of ISpruitl £l998) predicted 
very efficient angular momentum transport through the 
magnetic field coupling of core and outer envelope, the 
new quantitative formulation of magnetic field generation 
and magneti c field induce d angular momentum transport 
developed in ISm-uitl ll2002h leads to a more limited angular 
momentu m transport e f ficienc y. First calcul ations for mas- 
sive stars l|Heger et alll2005|) including the lSpruitl (|2002^ 
magnetic field physics show these stars evolving far from 
rigid rotation. Similar results are expected for low mass 
stars (Langer 2005, priv. comm.), which would likely bring 
the predictions for the white dwarf rotation velocities in 
agreement with the new observations. 

Our method could in principle be improved further, 
if that is necessary to distinguish between theoretical as- 
sumptions. The main reason for the current uncertainty 
is the unknown effective resolution due to rather large 
slits used in the spectroscopy. With a narrower slit and 
well controlled resolution an accuracy of 1-2 km/s could 
be reached, but the observations would need a significant 
amount of observing time on a very large telescope. 
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Table 2. Rotational velocities. When several spectra of 
the same object are of comparable quality, the results 
do not differ significantly and we therefore give only the 
result from the best spectrum. The only exception is 
WD2326 + 049, a known ZZ Ceti variable. [Ca/H] is the 
Ca abundance as determined from fitting the line pro- 
file, which is sometimes slightly different from that in 
the original source (see Table ^ the corrected EW for 
HS 0047+1903 used for the abundance in this table is 77 
mA) . V2 — v sin i is the rotation velocity determined using 
either the nominal resolution of the spectrographs (34000 
for Keck, 18500 for SPY), or a slightly higher resolution if 
the line depth cannot be reproduced at nominal resolution. 
In the latter case the values for the resolution are given in 
column Res. Errors are 1 a errors. V\ = (v sin i) max is an 
upper limit using the optimum resolution as explained in 
the text 
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Fig. 3. Seven objects fitted with Lorentz profiles. See text 



8 



L. Berger et al.: Rotation velocities of white dwarfs 



Table 3. Upper limit to the rotational velocities for the 
objects with deeper observed line profiles than is possible 
with a non-rotating model. The limits were determined 
using a Lorentz line profile with central intensity zero and 
equivalent width (EW, in niA) equal to the observed spec- 
trum 
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